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I. BACKGROUND 


A. INTRODUCTION: STATEMENT OF PROBLEM 

The rapid advancements in the field of miniaturization 
of electronic components is well-known. Technology has 
brought us from the large-scale integration (LSI) device to 
the very large-scale integration (VLSI) device in just a few 
short years. Since the first 1 kilobit RAM semiconductor 
device was introduced in 1970 the number of memory cells on 
a single device has grown to 64K while the size of the device 
Ma@eedecreased. One such device, the 64K RAM Micron chip pro- 
mecea by Tiny Micron Technology, Inc. 1s only 22,000 square 


ms, Or approximately the size of one of these printed 


letters [1]. Yet microminiaturization advancements are far 
from over -- 256K and greater RAM chips are on the immediate 
mor 1Z0n. 


As in most rapidly developing technologies there are var- 
10us problems associated with the advancements that have been 
Seagesin the Miniaturization of electronic devices. One of 
mie primary problems that has evolved is that of suitably 
cooling components so that they will perform reliably over 
TOng periods of time. For though the size of the individual 
component has shrunk considerably its junction temperature 
has remained fairly constant. It follows that the decrease 


in device size has been kin to an increase in heat flux. 


eZ 


To get an idea of the magnitude of this decreased size/ 
increased heat flux relationship one only has to heed the 
Meras Of Chu [2]. He asserts that the heat flux of a micro- 
@ercuit chip 1S but two orders of magnitude less than the 
heat flux of the sun ale vace ug! Watts per meter squared). 
Further, the microcircuit chip must dissipate its energy at 
strict maximum junction temperatures of between 125 to 150 
degrees C. The integration of a multitude of devices on a 
memgie printed circuit board leads to an increase in the 
temperatures in which the device operates. The high heat 
mex Of the microcircuit chip coupled with the increase in 
environmental temperature very often leads to the problem 
of component failure resulting from inadequate cooling. It 
is this problem that has recently prompted extensive research 


into the techniques of cooling microelectronic devices. 


Pee DESIGN CONSIDERATIONS: PROPOSED SOLUTIONS 

Many papers have been published outlining alternative 
methods of thermally controlling microelectronic devices. 
Hannemann [3] describes the thermal considerations that must 
be taken into account at each level (subassembly, assembly, 
and system levels) of the device design process. Kraus, et 
al [4] assert the need for thermal analysis even after 
design. They propose a strategy that yields conservative 
estimates of the heat flux associated with electronic cir- 


Sluitry. Further, Bar-Cohen, et al [5] indicate that numerical 





analysis of the heat transfer problem should be relied upon 
to cope with the analytically intractable nature of even 
standard packaging configurations. 

Once the initial design step has been completed a reli- 
able cooling scheme must be employed. A number of such 
schemes including conduction, forced air cooling, and direct 
liquid cooling have been proposed by many authors. Bergles, 
et al {6] contend that of the many techniques available, 
direct liquid cooling may be the most desirable mode of heat 
transfer when considerations such as weight, size, or cost 
are taken equally into account with the heat transfer coef- 
ficient. Simons and Moran [7] support this view by stating 
moat liquid coolants offer the capability to achieve heat 
transfer coefficients one to two orders of magnitude higher 
Phan that of air. 

Many authors have offered alternative methods of utilicz- 
ing immersion cooling. One such scheme is proposed by Kraus 
and Bar-Cohen [8] who analyze submerged condensers. They 
describe a model whereby the heat from the immersed elec- 
tronic devices is removed via bubble pumped augmented con- 
vection interacting at the submerged condenser surface. 

iimtseseldy £ocuses Om direct liquid cooling by natural 
convection. More specifically, since most microelectronic 
devices are seated on printed circuit boards packaged very 
near each other, investigation is limited to natural convec- 


tion in a small liquid filled rectangular enclosure. One 
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Mevantage Of Natural convection cooling is that there are no 
pumps, pipes, Or any other equipment within the enclosure; 


the manufacture of a relatively small module may be assured. 


C. IMMERSION COOLING: ANALYTICAL STUDIES 

Numerous analytical studies of natural convection in rec- 
tangular enclosures have been made since 1966 when Wilkes 
ena Churchill [9] initiated finite difference analysis for 
such configurations. This work was followed by many others 
Micluding Szekely and Todd {10] who, in 1970, presented a 
finite difference analysis of natural convection in a rec- 
tangular cavity in which two opposing vertical walls provide 
the temperature gradient. Comparing their numerical results 
With experiment they concluded that within the Rayleigh num- 
ber range of sixteen hundred to sixty-thousand the thermal 
response of a system is approximately independent of the 
Rayleigh number. They further noted a trend for secondary 
cells to develop once the Rayleigh number exceeded three 
hundred thousand. 

i 1975 Chu, et al [11] studied, via finite difference 
methods and experiment, the natural convection in a rectan- 
gular enclosure in which the heat source was a long, hori- 
zontal strip located in one vertical wall. Focusing on 
heater location, heater size, and enclosure size a conclu- 
Sion was reached that as the Rayleigh number increases the 


heater location for maximum circulation shifts downward but 
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that heater size seems to have little effect on circulation. 
They contend that secondary cells develop at the upper sur- 
face as the height of the enclosure is increased. 

Markatos and Pericleous [12] have recently published a 
study of buoyancy-driven flow in a square cavity with dif- 
ferentially heated sides. Secondary cells were also observed 


Metneir study. 


D. IMMERSION COOLING: EXPERIMENTAL STUDIES 

Megerlin and Vingerhoet [13] collected natural convection 
data using fluorocarbon dielectric fluids in an enclosure 
mee@eled after a microelectronic medule. One significant con- 
clusion reached was that cold plate orientation has a signi- 
ficant effect on the heat transfer -- a fifty percent increase 
Bemautained using a horizontal vice vertical cold plate. An 
experimental investigation using an undivided and partially 
divided water filled two-dimensional rectangular enclosure 
was performed by Nansteel and Greif [14]. 

In 1972 Baker [15] published the results of his experi- 
mental work in which he used resistors immersed in Freon. 
He concluded that the heat dissipation obtained with Freon 
1s at least four times greater than that obtained using the 
Same assembly with stagnant air. Lin and Akins [16] investi- 
gated natural convection heat transfer in a fluid filled cube 
where all six wall temperatures were subjected to a step 
change in temperature. Their conclusion was that a variety of 


circulation patterns evolve depending on the size of the cube. 
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In 1980 Ostrach, et at [17] published the results of an 
experiment that used small particles introduced in silicone 
Oils that filled various sized rectangular enclosures with 
Mmetght to length Patios less than unity. Secondary cells 
were evident in their study. Finally, an experiment was per- 
formed by Bohn and Kirkpatrick [18] which studied high Rayleigh 
number natural convection in a cublical enclosure filled with 
deionized water using various combinations of heated walls. 
They concluded that heating from below along with Penmikeaneous 
vertical wall heating and cooling enhances heat transfer. 

The above discussion is not intended to be a definitive 
Symoosis Of all works in the design considerations of micro- 
electronic devices, nor in the numerical analysis and experi- 
ments conducted on natural convection in enclosures. References 
Cited are merely intended to be representative of the works 
that have been done in these areas. They serve to provide an 
Seceilent point of departure for the remainder of this study. 
For a more extensive review of the works done concerning natu- 
ral convection in rectangular enclosures the reader is referred 


morcatton [19]. 


mee OBJECTIVES 

The objectives of this study are: 

1. to design and build an apparatus that approximates a 
liquid filled microelectronic package with appropriate printed 


circuit board spacing that is cooled by heat exchangers acting 


if 





as the top and bottom horizontal surfaces. The heat source 
Meco be a rectangular protrusion modeled geometrically after 
a 40-pin dual-in-line package (DIP) that is attached to and 
horizontally traverses the length of one vertical wall. 

Z2. to observe the steady state natural convection flow 
moeche apparatus using a variety of heater and heat sink 
memperatures as well as various locations of the heating 
erement . 

3. to analyze data obtained and determine Nusselt and 


Rayleigh numbers for each apparatus configuration. 


This study is intended to serve as a basis for future 
meat transfer experiments in natural convection using vari- 


Sus enciosure sizes and multiple heated protrusions. 
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A. GENERAL DESIGN CONSIDERATIONS 

The flow visualization apparatus (Fig. 1) consists of a 
small liquid filled rectangular enclosure with a heater 
modeled geometrically after a 40-pin DIP protruding into 
the enclosure from one of the vertical walls. This con- 
Preuration is meant to simulate a microelectronic component 
See printed circuit board. The top and bottom horizontal 
surfaces of the enclosure are such that their temperatures 
can be varied independently. Heat exchangers (Fig. 2) used 
an conjunction with a circulating bath act as these horizon- 
tal surface. A cartridge heater housed in a stainless steel 
cover (Fig. 3) acts as the 40-pin DIP model. This assembly 
is hereafter referred to as the heater. A calibrated resis- 
tor put in series with the heater and its power supply allows 
meee input power to be accurately determined. Eight thermo- 
memples are placed throughout the enclosure for temperature 
measurements. Two are placed on each heat exchanger and 
four on the cartridge heater cover. 

Several criteria were established to be used as guidelines 
for design. The following is a list of these criteria and 


maeir implications. 
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1. Dimensions 
As previously noted, the heater was to model geomet- 
Prcally a 40-pin DIP. Micrometer measurements of the height 
Mad width of a packaged 40-pin DIP were made and the rectan- 
gular heater cover is made to approximate these dimensions. 
Further, the enclosure's height was to be 6 times the verti- 
cal dimension of the heater cover and its width was to be 
between 3 and 4 times the heater cover's horizontal dimen- 
sion. This provides a good approximation to the spacing of 
a vertical printed circuit board placed in an electronic 
module. 
@. Heater Location 
The apparatus was to be made such that the heater 
could be raised or lowered on the vertical wall and that the 
width of the enclosure could be made to vary from 1 to 2 inches. 
This would provide a wide range of enclosure sizes and heater 
locations for investigation. It was determined that the 
ability to move the vertical wall on which the heater is 
attached is better for preserving fluid-tight integrity and 
smoothness. Machine milled interconnecting parts are success- 
fully used to meet this and the variable chamber width 
meqguirements. 
ee visud lization Technique 
The thymol-blue pH indicator flow visualization tech- 
nique proposed by Baker [20] was to be used. This is a non- 


intrusive technique that uses a pH indicator in distilled 





water. When a voltage is applied across the fluid it percep- 
moply darkens near the negative electrode due to a local 
Eqange in pH induced by proton transfer. If there is fluid 
motion away from the cathode, as there is in this experiment, 
the darkened fluid traces the path of motion before dissipat- 
mig Dack to the original lighter color and neutral pH of the 
Bulk fluid. In addition to its non-intrusive quality no 
density differences occur within the fluid during voltage 
meolication since the pH indicator remains an ion in solu- 
muomen the specifics of the fluid preparation are discussed 
in Chapter III. After much testing it was determined that 
the use of this technique requires any metal coming into 
contact with the liquid solution (e.g., heat exchangers, 
heater) to be relatively inactive with respect to the 
Galvanic series. This prevents corrosion as indicated by 
O'Conner [21]. To this end, a stainless steel shim strip 
running the length of the vertical wall opposite the heater 
is used as the cathode. The heater acts as the anode. 
This arrangement provides excellent visualization of the 
two dimensional flow of the darkened solution. 
4. Other Considerations 

One eighth inch plexiglas is used as the vertical 
Viewing wall to minimize optical distortion. Further, all 
vertical walls are heavily wrapped with foam rubber insula- 
tion to minimize heat loss via conduction through them 


during experimental operation. 
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B. COMPONENTS 
Peete d tel 

The heater assembly (Fig. 3) consists of an 8 inch 
momen 1/2 inch wide, i inch high, and 3/16 inch thick 304 
Stainless steel channel that houses a narrow rectangular 
cartridge heater backed by phenolic insulation that runs 
the length of the channel. The cover has four 1/16 inch 
deep slots located on the horizontal faces 2.75 inches from 
the ends of the channel. In each of these slots a copper- 
Genstantan thermocouple 1s spot welded nearest the cover's 
vertical face. The thermocouple wires run out the back of 
Biemsiots into the channel. The slots are filled with steel 
epoxy. Three evenly spaced screw holes are tapped into the 
top. and bottom of the cover so that the entire assembly can 
be attached to a vertical wall in the enclosure. Priot 360 
installation the cover is polished to a mirror-like finish 
using 180 to 600 grit special silicon carbide grinding paper 
and then cleaned with common rubbing alcohol. 

2. Heat Exchangers 

The two aluminum heat exchangers measure 1.5 by 2.5 
by 11.0 inches. The top heat exchanger is seen in Fig. 2. 
They each have three rectangular channels separated by 1/8 
inch walls. The channels are machine milled into a 1.375 by 
2.5 by 9.25 inches solid piece of stock aluminum and then 
covered with a 1/8 inch thick plate secured by ten no. 4 


screws. Slots 1/8 inches deep are milled 8.0 inches apart 
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On the side opposite the cover plate to accept vertical walls. 
One slot is 3/16 inches in width and the other 3/4 inches in 
width. Both slots traverse the width of the heat exchanger. 
In addition, two 1/16 inch deep slots are milled 2.75 inches 
from the outer slots and they run to within 1 inch of the 
other side of the exchanger. These last two slots have 
copper-constantan thermocouples spot welded into them and 
are covered with steel epoxy. A1.5 by 8.0 inch 0.001 gage 
stainless steel shim facing is epoxied in between the two 
outer slots on the side of the exchanger opposite the cover 
Beate. it 15 positioned flush against the corner of the 
Meat Exchanger that is one inch from the thermocouple slots. 
Mmemcoetacing 1s the relatively inactive portion of the ex- 
changer that acts as the horizontal wall of the chamber. 

The 1.5 by 2.5 by 0.875 inch headers are also machined 
from single pieces of aluminum stock. They each have a 1/4 
mmemeod, 1/8 inch 1d exterior pipe for coupling with a cir- 
culating bath and a 2 square inch hollow interior positioned 
SO as to get maximum flow in each exchanger channel. The 
headers are held to the main body of the exchangers by a 
Single no. 6 screw. All seams are fitted with epoxy prior 
to assembly. 

Both heat exchangers were tested for temperature 
distribution and leakage prior to use in the experiment. 
They were coupled to a circulating bath and temperature 


measurements were taken using the two installed thermocouples 
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and an additional similar thermocouple that was randomly 
moved about the exchanger surface. The thermocouples were 
run to a digital pyrometer with an accuracy of 1 degree 
Fahrenheit. It was determined from this test that the heat 
exchangers are, upon reaching steady state, at a constant 
temperature throughout; i.e., the pyrometer reading was the 
same for each thermocouple. 
3. Test Chamber 

The test chamber seen in Fig. 1 forms a rectangular 
enclosure measuring 1.5 inches in width, 8.0 inches in depth, 
and 6.0 inches in height. The top and bottom horizontal sur- 
faces of the enclosure consist of the stainless steel shim 
facing of the heat exchangers. Three of the vertical walls 
are entirely plexiglas. The wall holding the heater is 10.5 
anches wide, 24 inches high, and 11/16 inches thick. The 
heater is secured to the center of this wall with six no. 4 
screws and rubber o-rings that are counter sunk into the 
Poin «6oln addition, a 1/2 inch diameter hole is located 1 
inch from one end of the wall to accommodate the cartridge 
heater power leads and thermocouple wires from the heater 
cover. The wall has two 1/4 inch deep slots, one 3/16 
inches in width and the other 3/4 inches in width, 8 inches 
apart and 14 inches long running parallel along its height. 
These slots receive the adjacent vertical walls. They also 
allow the wall and subsequently the heater to be raised or 


lowered while maintaining fluid tight integrity. 
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Memeonner two plexiglas vertical walls are the ends 
mac comprise the 1.5 by 6.0 inch dimensions of the enclosure. 
They are both 2 inches by 6.25 inches but one is 11/16 inches 
thick while the other is 1/8 inches thick. The thinner wall 
is used for visualization and photography and was chosen to 
be thin in order to minimize optical distortion. Both walls 
slide into the 1/4 inch slots of the two larger adjacent walls 
and into the 1/8 inch deep slots on the heat exchangers. 

The final vertical wall of the enclosure is a 10.0 
gnch by 6.0 inch plexiglas frame that house a glass plate. 

The plate is meet ied with epoxy and will be used for future 
feuadies in conjunction with laser doppler velocimetry. Like 
the other large wall it has two 1/8 inch deep slots that 
receive the adjacent walls. The wall has lips that slide 
into the heat exchanger slots as well as two 1/16 inch id 
Standpipes made of 1/8 inch stainless steel tubing. One 
Standpipe 1s located near an upper corner and the other is 
diagonally opposite. They are used for filling/emptying the 
enclosure and relieving internal pressure during experimen- 


@at Operation. 


C. ASSEMBLY 

Once the heater has been attached to the large wall the 
apparatus 1s assembled with the exception of the glass wall 
framed by plexiglas. Dow Corning 732 RTV adhesive/sealant 
1s applied externally to all seams not involved with visuali- 


Zation, to the power lead/thermocouple hole, and to the heads 
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Sireall heater attaching screws. A thin coat of Dow Corning 
3140 RTV coating and silastic is used on all internal seams 
and on all of the external seams of the smaller walls. The 
glass wall framed by plexiglas is then fitted into place 
along with a taut 1/4 inch wide strip of 0.001 gage 302 
Stainless steel shim. The strip is centered internally on 
the wall and runs its height. A few inches of the strip 
traverse the seams between the wall and the heat exchangers 
(both top and bottom) to the outside of the test chamber. 
This strip acts as the cathode for the thymol-blue visuali- 
zation technique. The sealants are then applied to all 
remaining seams and the apparatus 1s tightened by specially 
made clamps. Fig. 4 is a photograph of the assembled apparatus. 
The vertical wall of the apparatus that holds the heater 
has to be removed, placed at the desired heater height, and 
then resealed between test runs. If removing the wall between 
Mest runs iS mot desired the apparatus can be turned upside- 


down to attain a different heater height. 


D. INSTRUMENTATION 
1. Power to the Heater 
[iewseater 1S rum in series with a precision resistor 
that was measured to have a 2.0317 ohms resistance on a 
Rosemount commutating galvonometer bridge. This arrangement 
allows a voltage measurement across the resistor to be made 


and, using its known resistance, the current through the 
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Figure 4. Photograph of Experimental Apparatus 
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circuit can be calculated. A measurement of the voltage 
across the heater is also taken. The product of the current 
and the voltage across the heater gives input power to the 
heater. Both voltages are measured using a Dana digital 
multimeter. 

Because of the 180 watt, 120 volt rating of the car- 
tridge heater two different power supplies are used in order 
to attain certain voltage levels. A Lambda regulated DC 
power supply is used for voltages up to 60 volts. A Variac 
voltage AC transformer is used for voltages greater than 60 
volts. <A concern about using the latter was the fluctuating 
voltage provided by house current and its possible effect on 
measurements. As it happened, temperature measurements re- 
mained stable. However, voltage measurements across the 
heater fluctuated up to 3.0%. Subsequently, calculations 
of input power to the heater are less accurate with this 
mode of operation than the power calculations made when the 
Lambda regulator is used. A better method when using the 
Variac would have been to put it in series with a line volt- 
age regulator in order to eliminate voltage fluctuations 
Caused by house current. 

ee voltage Across the Fluid 

The voltage source for the electrochemical technique 
is also a Lambda regulated DC power supply. Voltages across 
eac electrodes range from 3.0 to 8.0 volts depending upon 


the contrast of the fluid required for photography. The 
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voltage is measured with the Dana digital multimeter. The 
Beainless steel shim strip running the height of the verti- 
cal glass wall framed by plexiglas acts as the cathode while 
the stainless steel heater cover acts as the anode. Elec- 
trical access to the heater cover 1S via a wire spot welded 
to the bottom center stainless steel attaching screw that 
holds the heater assembly to the vertical wall opposite the 
cathode. 

Originally the heater was utilized as the anode and 
the top heat exchanger acted as the cathode. This arrange- 
ment proved unacceptable because flow paths could not easily 
be discerned from the viewing perspective. The inclusion of 
the shim strip as the cathode solved this problem by limit- 
ing the darkened flow to a single path located in the center 
Srethe enclosure. 

3. Temperature Measurement 

Eight, 30 gage copper-constantan (teflon/teflon in- 
Sulation) thermocouples dispersed throughout the enclosure 
are used for temperature measurements. Two thermocouples 
are located on each horizontal face of the heater cover and 
two on each heat exchanger. All thermocouples were fabricated 
using a Dynatech welder and calibrated on a Rosemount commu- 
tating bridge in conjunction with a platinum resistance tem- 


perature standard and calibration bath. 


The thermocouples run to a Newport digital pyrometer 


model 267 where temperature readings are made. The pyrometer 
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muves Teaaings to the meéarest degree Fahrenheit. If, during 
the measuring process, a reading fluctuates between two 
temperature extremes it is assumed to be the average of 

the two. Finally, the ambient temperature is measured with- 
in two feet of the apparatus with another thermocouple run 


mentne pyrometer. 
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Tiere CROMEN TA PROGEDURE 


A. ELECTROCHEMICAL FLUID PREPARATION 

Ms necea Carlier, the thymol-blue pH indicator flow vis- 
ualization technique as described by Baker [20] is used. 
Further, O'Conner suggests using thymol-blue sodium salt 
Vice pure thymol-blue (thymolsulphonephthalein) crystals in 
the solution to avoid filtering of undissolved crystals [21]. 
This advice is followed in the preparation of the fluid. 
First, 0.2 grams of thymol-blue sodium salt are dissolved 
in 1.0 liter of distilled water. This gives a 0.02% by 
weight light colored solution with no undissolved solids. 
Next, 0.1 normal sodium hydroxide is added drop by drop to 
mae Solution until it titrates to a dark blue color. The 
@ark color is an indication of a slightly basic pH. Finally, 
Maem sOlLULION 1S titrated back to its original lighter color 
by adding drop by drop 0.1 normal hydrochloric acid. The 
fluid is now slightly acidic and very close to its neutral 
pH value of 7. After preparation, a small amount of fluid 
1s poured into a shallow glass container. Wires attached to 
al.5 to 7.5 volt dry cell battery are used as the electrodes 
and the fluid is tested for local color change. All tests 


indicate adequate color change. 
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B. APPARATUS PREPARATION 

After the sealants are cured, 1/4 inch tygon tubing is 
attached to the heat exchangers using stainless steel hose 
Sramps. The exchangers are then interconnected in series to 
merei:rigerated circulating bath placed six feet from the 
M@pparatus. the top heat exchanger is placed upstream of the 
bottom exchanger. Common tap water is used as the refriger- 
ated circulating coolant. The apparatus is then laid hori- 
zontally on its largest wall. 1.5 feet of 1/8 inch tygon 
tubing is attached to each standpipe and the enclosure is 
ilied with the thymol-blue solution via a simple siphon 
System. Manual manipulation of the apparatus and standpipe 
tubing allows most bubbles to be driven out of the enclosure. 

A specially made steel leveling stand with four height 
adjustable legs and a slot for the largest vertical wall had 
been made to hold the apparatus. The stand is placed on a 
piece of plywood under which lies two inches of foam rubber 
insulation. This system serves to damp out any vibrations 
Occurring in the laboratory. The apparatus is placed verti- 
cally upright onto the leveling stand and secured in place. 
The two larger vertical walls are then layered with 3/4 inch 
thick insulation. After insuring that the view through the 
test chamber to a white background is unobstructed, 3/4 inch 
thick insulation is placed on the remaining vertical walls. 
The apparatus is then leveled both horizontally and vertically 


by means of a bubble level placed on the top heat exchanger 
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in two orthogonal directions and by concurrent adjustment of 
the stand legs. This leveling process is crucial for ensur- 
ing horizontal and vertical alignment of the chamber. Fig. 35 


is a photograph of the assembled apparatus and instrumentation. 


fee LES! PROCEDURE 
1. Heater Location 

Test runs were made for three heater locations: (1) 
the top horizontal face of the heater located 1.0 inch from 
the top of the enclosure, (2) the heater centered on the 
vertical wall of the enclosure, and (3) the bottom horizontal 
face of the heater located 1.0 inch from the bottom of the 
emelLOSure. 

a) inl tent nMmenmc Settings 

Temperature measurements of all thermocouples are 
made to insure that the apparatus and its liquid contents 
Pema 4a Singleetemperature. ff this criteria is met the 
circulating bath to the heat exchangers is turned on to a 
Peeadetermined temperature. Shortly thereafter the power 
Supply to the heater is turned on and the voltage slowly 
increased to another predetermined value that results in 
the approximate desired temperature difference between the 
exchangers and the heater. The values of the settings for 
the circulating bath temperature and the voltage to the 
heater are determined soley from experience garnered from 


a multitude of test runs. 
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S inserumemt Readings 


Thermocouple temperature measurements and voltage 
measurements across the heater and resistor are taken every 
twenty minutes until successive temperature measurements 
memain unchanged to within the accuracy of the pyrometer. 
Measurements are then taken every ten minutes. It is assumed 
that steady state is attained within the enclosure if no 
thermocouple measurement changes by more than one degree 
during any ten minute period. This process takes between 
maomto four hours depending on the initial temperature dif- 
ference between the heater and the exchangers. The smaller 
temperature difference takes the least amount of time. 

Once steady state is attained the power supply that 
creates the voltage across the electrochemical fluid is 
turned on. Using a digital multimeter as an indicator, the 
voltage is increased from 0.0 to 3.0 volts. The insulation 
1s then removed from the two smaller vertical walls and the 
fluid voltage is varied until a visually suitable path of 
darkened fluid ensues. 

4, Photographic Techniques 

A Nikon Ftn with a 50 mm f2 lens is used for all black 
and white photography. A Canon Ftb with a 50 mm f1.4 lens 
1s used for all color photography. Black and white photo- 
graphs are taken with Kodak Tri-x pan asa 400 film; color 
photographs are taken with Kodak Ektachrome asa 64 daylight 


film. Black and white photographs are taken first. Both 
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Cameras are secured to a tripod. The Nikon is placed between 
meandeiS inches from the 1/8 inch thick plexiglas vertical 
wall. A combination of f-stop settings of 8 and 11 coupled 
with shutter speeds of 1/60 and 1/125 of a second are used. 
The Canon is placed at a distance of between 12 and 15 inches 
from the viewing wall. Combinations of f-stop 2 and 4 with 
shutter speeds of 1/4 and 1/15 of a second are used with 

eae Canon. 

Two L-shaped pieces of cardboard are attached together 
feeds tO create a rectangular opening approximately the size 
of the small vertical wall. This is placed between the camera 
and the enclosure. A cloth is hung from this cardboard aper- 
ture and draped over the camera during photography to elimi- 
Nate camera reflection on the plexiglas. Another larger 
piece Of cardboard covered with white paper is placed 1 foot 
from the rear small vertical wall to act as a backdrop. 
Finally a 650 watt colortran multi-6 set lamp is positioned 
2 feet from the white cardboard backdrop and 1.5 feet from 
the apparatus. The lamp is directed towards the white back- 
drop and turned on only when photographing the flow. 

A major concern was that the heat generated by the 
lamp would disturb the flow. However, due to the insulation 
remaining between the lamp and the enclosure and the short 
interval of time the lamp is on the viewed flow remains two 
dimensional. This scheme provides an excellent view down the 
length of the enclosure to a white background that adequately 
contrasts the darkened fluid. 
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D. DATA ANALYSIS 
There are six length dimensions within the enclosure, any 
one of which could be chosen as the characteristic length, 
L, for determining the Nusselt number and Rayleigh numbers. 
The 1.0 inch vertical face of the heater is chosen since it 
is the surface that not only adds the majority of energy to 
the system but it also is oriented in the direction of the 
gravitational driving force of free convection. The proper- 
ties of the fluid are evaluated at the average of the hot 
surface (heater) and cold surface (heat exchangers) temperatures. 
1. Determination of Nusselt Number 


The Nusselt number is defined as follows: 


nae 

Nu = aes (ay) 
Memevaluate the heat transfer coefficient an estimate has to 
be made of the heat lost through the vertical walls of the 
enclosure. <A parallel thermal resistance circuit for con- 
duction through the three vertical walls in contact with the 
heater is employed for this purpose. Heat loss through the 
vertical wall opposite the heater, all contact resistances, 
and conduction via the heater assembly attaching screws are 
assumed negligible. Knowing the total thermal resistance 
meOomeconduction the heat loss is calculated using a variation 
Of Fourier's Law of Conduction: 


We = ue 


Q = =.= 
cond ee 


(Za) 
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where 
b 


Reond ~ EK (2b) 
iiestength of the conduction path, L, for all materials is 
known. The thermal conductivity, k, for each material is 
determined using a table of properties [22]. The area, A, 
1s assumed to be that of the cartridge heater pe LpemidTeurar 
mMomene COnduction path. Finally, the temperature difference 
is assumed to be the difference between the ambient, Ue 


pide tne heater cover temperature, T Gis Lurehen meglects 


s? 
the convective thermal resistance on the outside of the in- 


sulation. In the event ee 1s less than T heat loss is 


4? 
assumed to be zero. 

After calculating the energy lost via conduction and 
having already determined the energy into the system (Q.) 


from the voltage measurements across the heater and the pre- 


cision resistor a simple energy balance is used to determine 


the energy convected into the fluid Care): 
Qconv * Qin - Rcona (3a) 
where 
V, V> 
in ep) 
and 


< 
i 


1 voltage across heater 


V, = voltage across resistor 
20527 Ohms 


ye) 
il 
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Then, using Newton's Law of Cooling, 


Qconv = (T. T,) (4) 
h is calculated. Ty, is chosen to be the temperature of the 
heat exchangers which are the ultimate heat sink for the 
enclosure. The exterior surface area of the heater cover 
is A. With all variables in equation (1) now known the 
Nusselt number is calculated. 

To ensure the accuracy of this method of analysis a single 
iteration is performed on the Nusselt number by use of another 
thermal resistance for conduction Circuit. This time a single 
circuit with a lone resistance is run through the steel heater 
cover to the cartridge heater with the intent of determining 
meewecartridge heater temperature. Using the heater cover 
temperature Cla) 5 length of the conduction path through the 
heater cover (L), the thermal conductivity for 304 stainless 
steel (k..), the exterior area of the heater cover (A), and 
the previously calculated convection energy into the enclo- 


sure (Q ), 1t is a simple matter to determine the cartridge 


COnV 


heater temperature (T)) from a rearrangement of Fourier's Law 


S@eeconduction: 


LQ 
— conv , + (5) 


Ss ; 
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i in equation (2a) is then set equal to Tj; and a single 
meeration 1s performed to re-evaluate the Nusselt number. 
Iterative analysis of Nusselt numbers for all test runs 
mm@ove to fall within 3.03 of the initial calculation thus 
eliminating the necessity for further iterations. A sample 
calculation is contained in Appendix A. 

fee Levermination of Rayleigh Number 


One definition of the Rayleigh number is as follows: 


Ra = (Gr) (Pr) (6a) 
where 
3 
es I Cl Se 
aa ee (6b) 
v 
Pr= = (6c) 


In the calculation of the Grashof number (Gr), the same 
characteristic length (L) is used as that used to determine 
the Nusselt number; that 1s, the 1.0 inch vertical face of 
the heater cover. In addition, the same temperature differ- 
ence is used as that in equation (4) (i.e., the average 
heater cover temperature minus the average heat exchanger 
temperature). The remaining terms in equation (6b) are eval- 
uated at the average temperature of all thermocouple measure- 
ments by a cubic polynomial that had been fit to data from 

a table of properties [23] for temperatures of water ranging 


from 50 to 100 degrees Fahrenheit. 
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The Prandtl number (Pr) is also evaluated by a cubic 
polynomial that had been fit to data from the same table 
of properties, at the same temperature, and for the same 
working fluid. The Rayleigh number is then calculated. 

A sample calculation is contained in Appendix A. A sample 


gee rtainty calculation 1s contained in Appendix B. 
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IV. RESULTS 


eee LSUAL 

Photographs of the natural convection flow within the 
enclosure for a representative number of heater locations 
and Rayleigh numbers are presented in Figs. 6, 7, and 8. 
meee 6 depicts flows with the top horizontal face of the 
Memeer located 170 inches from the top of the enclosure. 
Beeme 7 depicts flows with the heater centered in the en- 
mesure. Finally, Figure 8 depicts flows with the bottom 
horizontal face of the heater located 1.0 inches from the 
bottom of the enclosure. 

Some general observations can be made about the natural 
convection flow within the enclosure. First, the flow 1s 
Suar-celled. There is a buoyancy driven upper cell and a 
Shear driven lower cell. The location of the cell separa- 
tion layer is dependent upon the heater location within the 
ememosure; aS the heater 1s raised so, too, is the layer. 

Secondly, the flow of the buoyant driven upper cell fol- 
lows the geometry of the enclosure. The fluid comes off 
Miemvertical face of the heater, up the vertical wall to 
which the heater is attached, across the cold heat exchangers, 
then down the vertical wall opposite the heater. Finally, 
the fluid moves along the cell separation layer where it 


interacts viscously with the lower cell and back towards the 
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heater. No definitive comments can be made about the total 
motion of the fluid in the lower cell; only a very slow moving 
flow was observed slightly below the separation layer. This 
flow was crossing the width of the enclosure towards the ver- 
tical wall holding the heater. 

The photographs also show flow being generated into each 
meee | 6This is the result of a recurring problem encountered 
using the Baker flow visualization technique. Baker suggests 
Poplying a2 maximum of 5.0 volts across the fluid when elec- 
meode spacing is 0.5 cm. These limitations preclude genera- 
mpom ot bubbles near the cathode. Utilizing the limits imposed 
and a network of cathodes the ensuing flow can be seen as it 
moves away from the area of interest. In this experiment a 
fees Cathode is used. The area of interest is the entire 
pmeeosure. Hence, it 1s necessary to apply up to 8.0 volts 
in order to observe the flow throughout the entire enclosure. 
Bubbles are therefore generated that cause convective mass 
Mmeamoter up the stainless steel shim strip. The buoyant flow 
Circulation of the upper cell collides with the mass transfer 
associated with the convected bubbles resulting in the intru- 
meommot flow into the cells. Efforts were made to minimize 
Mimamcttrect by careful adjustment of the voltage during test 
mums . 

iigescdistinet ditference between the photographs of Fig. 7 
Slearly illustrates the advantage of utilizing a stainless 


Seee shim strip as the cathode for the visualization technique 


49 





As previously noted, the heater was originally utilized as the 
anode while a heat exchanger acted as the cathode. This pro- 
vided an excellent darkening of the fluid but also proved in- 
adequate in terms of flow visualization (right photograph of 
pees’). The inclusion of the shim strip as the cathode allows 
visualization of the two dimensional flow that occurs in the 


middle of the enclosure (left photograph of Fig. 7). 


B. QUANTITATIVE 

Table I summarizes the results of calculations using data 
obtained for each test run. Based on these calculations the 
Nusselt and Rayleigh numbers were plotted on a log-log graph 
(Fig. 9). Since the temperature measurement uncertainty in- 
herent in the experiment produces Nusselt and Rayleigh number 
uncertainties ranging from five to twenty percent (see Appendix 
meet is difficult to draw definitive conclusions from the 
data. However, there does appear to be a trend suggesting 
that as the heater is raised within the enclosure the Nusselt 


number decreases. 
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TABLE I 


Summary of Results 
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V. CONCLUSIONS 


The author has found no previous natural convection study 
faa ligquad tilled rectangular enclosure with a heated pro- 
trusion attached to a vertical wall, temperature variant 
horizontal walls, and the remaining walls adiabatic. However, 
some investigators using similar configurations have reached 


the same conclusions about to be presented. 


A. VISUAL 

Two cells occur within the enclosure during the natural 
convection flow. The upper cell is driven by buoyant forces 
and essentially follows the geometry of the enclosure. The 
Moen cell is shear driven. The location of the cell separa- 
tion layer is dependent upon the height of the heater in the 


enclosure. 


B. QUANTITATIVE 

Calculations indicate a trend for the Nusselt number to 
decrease as the height of the heater within the enclosure is 
mierecased. It is interesting to note that for the Rayleigh 
numbers calculated in this experiment (on the order of 10° 
to 70°) Chu, et al [11] plot this same tendency for an air 
filled square enclosure that has a heater imbedded in a ver- 


meeai wall and cooled horizontal walls. It is Chu's inves- 


tigation that most closely resembles this experiment. 
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VI. RECOMMENDATIONS 


A. IMPROVEMENT TO EXPERIMENT 
ie) npmanat Us 

Fluid leakage into the heater assembly occurred during 
the course of the experiment. This ee the cartridge heater 
Memtail. Replacement cf the cartridge heater necessitated dis- 
mantling the apparatus. To avoid the problem of cartridge 
heater malfunction, either a liquid immersible cartridge heater 
should be used or else an improved design of the heater cover 
Should be sought. 

though the phnetographs proved satisfactory it is fur- 
ther suggested that experimentation with various camera filters 
be done. The contrast of the observed flow may improve to the 
point where excessive voltage across the fluid may not be re- 
quired, thereby, eliminating the bubble generation effect. 

2. Instrumentation 

In order to lower the experimental uncertainty it is 
Suggested that a more accurate temperature measuring instru- 
ment be used for thermocouple readings. Further, to preclude 
the use of the Variac transformer a cartridge heater with the 
Same power (180 watts) but rated at 60 instead of 120 volts 
Should be used. This will allow all test runs to be made 
with the Lambda power supply as the power supply to the 


heater. 
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B. ADDITIONAL EXPERIMENTAL WORK 

It is suggested that the following areas of study be 
experimentally explored: 

1. Decrease the enclosure width to see if any additional 
cells develop within the enclosure. 

2. Observe the effects of independently varying the heat 
exchanger temperatures. In conjunction with this, do not 
mimlaw the Coolant temperature to go below fifty degrees Fah- 
renheit to avoid the density inversion that occurs in water 
at forty degrees. As can be seen in Table I, the heat ex- 
Saaneer temperatures for some test runs fell below this value. 

Pee cOmtinue £6 use the stainless steel shim strip as the 
cathode but seek methods of decreasing the effect of the 
bubbles generated when voltage is applied across the fluid. 
Alternatively, design a network of cathodes (as suggested by 


Baker) that encircles the enclosure. 


55 





APPENDIX A: SAMPLE CALCULATIONS 


A. DETERMINATION OF INPUT POWER 

itiresdimensions ©: the cartridge heater are 8.0 by 0.25 
Dy 0.725 inches. Table II is the data that was recorded 
for a representative run on 3 November 1983 with the heater 
located near the bottom of the enclosure. Refer to Fig. 1 


more thermocouple locations. 


TAB EE tal 


Test Run Data for 3 November 19853 


2 3 4 os 6 17 18 “8 ee 2 
40 40 40 40 73 77 #77 #=+%76 71 Sko) aod Lea 
Notes: (1) all temperatures in degrees Fahrenheit 


(2) all voltages in volts 


Using equation 3b the input power is calculated to be: 


ese 51) (1.42) 


fe foeet oe ht ne-Watt) = 139.54 Btu/hr 
(20sa7 wonims | 


B. NUSSELT NUMBER 
The heat loss via conduction through the vertical walls 


is determined by constructing a thermal resistance network 
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(Fig. 10). The thermal conductivity of each material 


(table I11) is determined from a table of properties [22]. 


RUBBER 
INSULATION PLEXIGLAS 


HEATER 
COVER 





CARTRIDGE 
HEATER 


PHENOLIC 


Figure 10. Thermal Resistance Network for Calculation of Losses 





(eS Ea 


Thermal Conductivity of Materials 


Material K (Btu/hr-ft-°F) 
rubber insulation O54 S 

plexiglas Os) 24 11 

phenolic OPOn 7s 

304 steel Seog 5 


Table IV is a summary of the individual resistances for Fig. 


10 that were determined using equation 2b. 


TABLE IV 


Thermal Resistances to Conduction 


Resistance Path(1) Material Value (°F-hr/Btu) 
Ry a plexiglas eos, 
R., a rubber insulation Zoe ol 
Rz D plexiglas 643.2 
Ry D rubber insulation BoO7 5 
Re Cc phenolic One 
Re e plexiglas PAVE 
Ro E rubber insulation S 0a 


Note: (1) conduction path through 
a. viewing wall 
b. wall opposite viewing 
c. wall to which heater is attached 
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The equivalent resistance for each path is: 


R, = 2679.6 OF-hr/Btu 
Rye o0er 2 OR-hr/Btu 
eee OF-hr/Btu 


Therefore, 


(R,) (Ry) (RQ) 
ond ee ee TOs °F-hr/Btu 
(R,) (RL) +(R,) (RL) +(R,) (RD 


Employing equation 2a, 


_ (75.75 - 40.0) 
(101.28) 


ota = 0.047 Btu/hr 


The Nusselt number is determined by first using equation 


5a to calculate the energy convected into the fluid: 


eeremn@lcd2 40>) 0-047) = 139-495 Btu/hr 


Then, the heat transfer coefficient is calculated as (equa- 


cmon 4 ) 


Ze 


Gee ele 4) 35 10 Btu/hr-£t+-°F 


a> -40.0) (16) 


ie = 


The thermal conductivity of water is determined from a table 


of properties [25], 


5S, 





ke = 0.3383 Btu/hr-ft-°F 


Finally, the Nusselt number can be calculated (equation 1), 


(S32 LO RICW/ 170) en 
(0.3383) 


Nu = 


An iterative check is performed by determining the car- 


tridge heater temperature (equation 5) as 


ie = (3/192) (139.495) Peis Se =877 96 -F 
(8.659) (1.33) 


Reapplication of equations 2a, 3a, 4, and 1 respectively with 


qT. set eCqual to Th yields 


mee o One 7 1.0) 


Q = = 0.0687 Btu/hr 
cond (101.285) 


Qeony = (139.54 - 0.0687) = 139.474 Btu/hr 


ae = JUS ie Cee = SS ees £0, Btoynr-tt 


(736-71 .0) (16) 


2 Of 


Poem the table of properties [23] the fluid conductivity 


remains 


ke = 0.3383 Btu/hr/ft/°F 
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Therefore, 


Nu = (35.10) (1/12) 8.64 
(OMS5.835)) 


The percent difference in Nusselt numbers after one itera- 


Mom iS Zero percent. 


See RAYLEIGH NUMBER 
Using the property table [23] at the average temperature 
of the heater and heat exchangers the Prandtl number and 


Grashof number can be calculated as: 


Pie = Goce 7 
Ae Sse 9) F/ecubic ft 
MY, 


PGs cen? 96 = 40.0) (1/1728) = 322.4 x 10° 


Therefore, the Rayleigh number can be determined using 


equation 6a: 


ne ea esi) = 2.68 x 10° 
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APPENDIX B: UNCERTAINTY ANALYSIS 


Uncertaintiesin the experiment are evaluated by employing 


the root mean square method of Kline and McClintock [24]. 


Table V is a summary of measurement and property values and 


their associated uncertainties: 


Variable 


gB/v- 


ie (2) 
V,(3) 
v4) 


Notes: 


(1) 
(2) 
(3) 
(4) 


TABINE Vv 


UniGertaimty in Variables 


Deere adie y: 
L258 


00 
5 
(aa 
720 
> 
0 


oe vo Qo oe 


0 
BOO) cin 
1/600 
0.0001 ohms 
05 
02 eevoOlEs 
yy Ole volts 
ORD TLV OES 


bes Lommumecen Eaunmit y 


assume last digit 


assume last digit 


[22] 
[22] 
[22] 
[22] 


resolution of measurement device 


assume last digit 


Ga oaeat 1 On 


resolution 
resolution 
resolution 


resolution 


all dimension measurements 


Oar 
Of 
Ot 
of 


on galvonometer 

measurement device 
measurement device 
measurement device 


measurement device 


all thermocouple temperature measurements 


Lambda regulator as power supply 


Variac as power supply 


G2 





The following variable uncertainties are determined using 
the aforementioned RMS equation on the data presented in 


Table II of Appendix A. 


2 2 2 1/2 
(wv,)(V,) 7 cwv,)(v,)_2 GRD CV) CV) 
ee ee 
2 2 2 V2 
(£0.01) (1.42), “, (0.02) (58.51), °, - (0.0001) (58.51) (1.42), *, 


O57 (2.0317) (2.0317)2 


0.983 Btu/hr 


, 2 ie 
WT. Z Ni oe 2 Ciena a. =" Te) 
eee) (-— ) + (——————_,-—_) ]} 
‘cond Saat nena 
2 2 2 1/2 
err) / 2) Cay 2) (0.00005) (2) 
ines) 8G) + i) | 
(101.28) (101.28) (101.28) 
= 6.017 Btu/hr 
2 21/2 
Se OnV 7 Lee ig CiiGlaee 
17 2 


eM 0983) - 1 


eo So ees Ctl, hr 
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eZ 


Z Z Z 
wu = (WBC) ” 4 WLI)” (Ke) Ch) LD)“ 


(ky) (k,) (k,)¢ 
2 2 ae 
(2.95) (0.083) )~, (0.001) (35.10))",- (0.001) (55.10) (0.083))", 
(0.338) (0.338) (Onsen - 
= 0.727 
; a 2 
eee (((WGr)(Pr))” + ((WPr)(Gr))°] 


2 172 


Wace 10, 08 322))- + (€0.0)(322.0 x 10°)) | 


memes x 107 


The above uncertainty calculations are intended to be 
mepresentative of the method used; calculations performed on 
every variable are not included. A summary of quantitative 
results for Nusselt and Rayleigh numbers and their respective 
Meer tainties in relation to heater temperature minus ex- 


Smameer temperature for all test runs is given in Table VI. 
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TABLE VI 


UiGepeatityeliuexperiment 


Wet Nu SOUNGemeainty Ra (x 10°) UNGER tainuy 
Heater Near Top 

el. 7 5 4.95 See 0 5 SS 17.40 
SUR aS Soul Jed | aes eee OS) 4 
Bo, 58 6.84 9S) g lols, Bas U0 
42.38 (eos U xyes Ou.J0:) SSS 
43.25 Teo8 ie 9 SrtOS Sy ite: 
ao. 00* hg ILS (59 oS ae: 
Heater Centered 

DHS a EMS Oe 2 MNRAS 
HO ZS 5.0.6 CASS es ile ON iS) Rous 
eee OO a0 ee. 2S 109 Issa eye: 
HES. 5.0 Gee 12S ess, 11.40 
mee! 3 eZ Des 8) 1.49 2) Re 
3) tes) U5 ons Se die 1.94 (Boat 
5) as) Ws eke: S97 oo Ce 
Bc. 00 8.46 Seles S00 6.59 
Be8,/5 8.46 ea be) 20 Deol 
So. 20 e268 Srey eee | O20 
40.75 SO or S208 Orca 
48.00* 8.46 SO Sel! eel /4 
open 50 * SS Las 5 OO STS 4.87 
Heater Near Bottom 

Ieee / 5 2.40 seas Be) rales 17.40 
iy. 00 6.66 Sp 5 1 1.40 Wars Sal 
ZS, 20 ees 9.43 lay ale Fee) 0) 
ao. / 5 8.64 8.41 2508 opto} /f 
Note: (1) * Indicates Variac used as power supply. 
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